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Abstract 

We consider the spherically symmetric, asymptotically flat 
Einstein- Vlasov system. We find explicit conditions on the initial data, 
with ADM mass M, such that the resulting spacetime has the following 
properties: there is a family of radially outgoing null geodesies where 
the area radius r along each geodesic is bounded by 2M, the time- 
like lines r — c G [0, 2M] are incomplete, and for r > 2M the metric 
converges asymptotically to the Schwarzschild metric with mass M. 
The initial data that we construct guarantee the formation of a black 
hole in the evolution. We also give examples of such initial data with 
the additional property that the solutions exist for all r > and all 
Schwarzschild time, i.e., we obtain global existence in Schwarzschild 
coordinates in situations where the initial data are not small. Some 
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of our results are also established for the Einstein equations coupled 
to a general matter model characterized by conditions on the matter 
quantities. 

1 Introduction 

One of the many striking predictions of General Relativity is the assertion 
that under appropriate conditions astrophysical objects like stars or galaxies 
undergo a gravitational collapse resulting in a spacetime singularity. This 
was first proven by Oppenheimer and Snyder [TO] who constructed a semi- 
explicit example of a homogeneous spherically symmetric ball of dust, i.e., 
of a pressure- less fluid, which under its self-consistent, general relativistic 
gravitational interaction collapses. During this collapse the scalar curvature 
of spacetime blows up at the centre of symmetry, and the geometry of space- 
time breaks down there. This is referred to as the formation of a spacetime 
singularity. An important feature of the Oppenheimer-Snyder solution is 
that during the collapse a two-dimensional spacelike sphere evolves which 
encloses the singularity and through which no causal curve, i.e., no light 
ray or particle trajectory, can pass outward. In this way the spacetime sin- 
gularity is isolated from the outside part of spacetime by a so-called event 
horizon, and the singularity cannot be seen or in any other way be experi- 
enced by observers outside the event horizon. This configuration was later 
termed a black hole. 

In the 1960s Penrose [20] proved that the formation of spacetime singu- 
larities from regular initial data is not restricted to spherically symmetric, 
especially constructed or isolated examples but is a genuine, stable feature 
of spacetimes. However, this result gives little information about the geo- 
metric structure of a spacetime with such a singularity. In particular, it is 
in general not known if every spacetime singularity which arises from the 
gravitational collapse of regular initial data is covered by an event horizon. 
Since the existence of so-called naked singularities (for which, by definition, 
the latter is not true) would violate predictability (it would not be possible 
to predict from the initial data what an observer would see if he could ob- 
serve a singularity), the cosmic censorship conjecture was formulated which 
demands that any singularity which arises from the gravitational collapse of 
generic regular initial data is indeed hidden behind an event horizon. The 
restriction to generic data means that naked singularities are allowed to oc- 
cur for a "null set" of the initial data. An important example where naked 
singularities do form for a null set, but for which cosmic censorship holds 
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true, is the spherically symmetric Einstein-scalar field system, cf. \10\ 111]. 
Actually the above is an informal statement of the so-called weak cosmic 
censorship conjecture |34t 12.1]; we will not be concerned with the strong 
version in the present paper. For a mathematical discussion and the defini- 
tion of the weak cosmic censorship conjecture we refer to |12]. 

To deal with this conjecture in full generality is out of reach of the 
present level of mathematics, but under the assumption of spherical sym- 
metry progress has been made in recent years. One important outcome of 
these investigations is that the answer is sensitive to which model is chosen 
to describe the matter. Christodoulou [7] showed that for dust, i.e., the 
matter model used by Oppenheimer and Snyder, cosmic censorship is vio- 
lated. On the other hand, in a series of papers Christodoulou investigated 
a massless scalar field as matter model and showed in 1999 that weak and 
strong cosmic censorship hold true for this matter model; see [llj and the 
references therein. 

In the present investigation the main example considered as a matter 
model is the so-called collisionless gas as described by the Vlasov equation. 
It is used extensively in astrophysics, cf. [6j, to describe galaxies or globular 
clusters which are viewed as large ensembles of mass points which interact 
only through the gravitational field that the ensemble creates collectively. 
In a relativistic context this leads to the Einstein- Vlasov system. All results 
available for this system support the following 

Conjecture: Weak cosmic censorship holds for the Einstein-Vlasov system. 

We mention explicitly that, in contrast to dust, small, spherically symmet- 
ric initial data launch global solutions, i.e., the solutions are geodesically 
complete and hence satisfy cosmic censorship, cf. |25j . Also, the numerical 
simulations [5], UHl [29] which treat large initial data support the hypothesis 
that naked singularities do not form in the evolution. We point out a further 
interesting feature of Vlasov matter observed in these numerical studies: In 
a one-parameter family of solutions which for large parameters, i.e., large 
amplitudes of the initial data, collapse to a black hole the smallest black hole 
always has a strictly positive ADM mass, i.e., there is a mass gap. For some 
other models, e.g. a scalar field, the mass of the black hole as a function of 
the parameter is continuous and arbitrarily small black holes can form, cf. 
)15] for a review. 

The aim of the present paper is to find explicit conditions on the initial 
data which ensure the formation of black holes. This class of initial data 
has the important property that, except for "boundary cases", properly re- 
stricted small perturbations of the data lead to solutions with the same 
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properties. In this sense the estabhshed behaviour of the solutions is stable 
and not restricted to especially constructed solutions or initial data, respec- 
tively. It turns out that some of our results can be formulated for a general 
matter model which satisfies certain specific assumptions, and in order to 
give a broader impact to our results we shall do so. At the same time we 
emphasize that the Vlasov matter model is the only one which is presently 
known to actually satisfy all the assumptions needed for our arguments to 
go through. 

As an interesting corallary to our main result we show that it is in fact 
possible to choose initial data for the Einstein- Vlasov system, which lead to 
formation of black holes, such that the solutions exist for all Schwarzschild 
time and all r > 0. We thus obtain global existence in Schwarzschild coordi- 
nates for initial data which are not small, and to the best of our knowledge 
this is the first global existence result in Schwarzschild coordinates for initial 
data which lead to gravitational collapse and formation of black holes. 

One aspect of our result is that there is a set of initial data which leads 
to gravitational collapse such that weak cosmic censorship holds. This point 
should be related to an earlier result by Rendall [32] , where it is shown that 
there exist initial data for the spherically symmetric Einstein- Vlasov system 
such that a trapped surface forms in the evolution. The occurrence of a 
trapped surface signals the formation of an event horizon. Indeed, Dafermos 
|13j has proved that if a spherically symmetric spacetime contains a trapped 
surface and the matter model satisfies certain hypotheses then weak cosmic 
censorship holds true. In [I4j it was then shown that Vlasov matter does 
satisfy the required hypotheses. Hence, by combining these results it follows 
that initial data exist which lead to gravitational collapse and for which weak 
cosmic censorship holds. However, the proof in [32] rests on a continuity 
argument, and it is not possible to tell whether or not a given initial data 
set will give rise to a black hole. This is in contrast to the explicit conditions 
on the initial data, together with the detailed asymptotic structure, that we 
obtain in the present work. In this regard it is natural to relate our results 
to those of Christodoulou on the spherically symmetric Einstein-scalar field 
system [8] and [9]. In [8] it is shown that if the final Bondi mass M is 
different from zero, the region exterior to the sphere r = 2M tends to the 
Schwarzschild metric with mass M. In Theorem 12.41 below we show that 
solutions of the spherically Einstein- Vlasov system, under certain conditions 
on the initial data, also converge to the Schwarzschild metric asymptotically. 
Furthermore, in [9] explicit conditions on the initial data are specified which 
guarantee the formation of trapped surfaces. This paper played a crucial 
role in Christodoulou's proof [11] of the weak and strong cosmic censorship 
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conjectures mentioned above. The conditions on the initial data in [9] allow 
the ratio of the Hawking mass and the area radius to cover the full range, 
i.e., 2m/r S (0,1), whereas our conditions always require 2m/r to be quite 
close to one. However, we believe that to understand gravitational collapse 
in the case of Vlasov matter the essential situation is when 2m /r is large. 
We thus hope that the results in the present paper will lead to progress 
on the general understanding of gravitational collapse and the weak cosmic 
censorship conjecture in the case of Vlasov matter. 

The Vlasov matter model has a further property to recommend it when 
compared to other matter models. For the Vlasov-Poisson system, which 
arises as the Newtonian limit of the Einstein- Vlasov system in a rigorous 
sense [261 [31], and which is used extensively in astrophysics, there is a global 
existence and uniqueness result for general, smooth initial data [17ll21j . This 
means in particular that any breakdown of a solution of the Einstein- Vlasov 
system can be expected to be a genuine, general relativistic effect such as a 
spacetime singularity and not only remainder of some bad behaviour which 
the matter model exhibits already on the Newtonian level. 

To be more specific, consider now a smooth spacetime manifold M 
equipped with a spacetime metric gap] Greek indices run from to 3. Then 
the Einstein equations read 

Gap = 8TTTa/3, (1-1) 

where Ga/3 is the Einstein tensor, a non-linear second order differential ex- 
pression in the metric ga/s, and Ta/s is the energy-momentum tensor given 
by the matter content (or other fields) of the spacetime. To obtain a closed 
system, the field equations (jl.ip have to be supplemented by 

evolution equation(s) for the matter (1-2) 

and 

the definition of Tap in terms of the matter and the metric. (1-3) 

It is often possible to specify conditions on ()1.2p and ()1.3p under which one 
can establish geometric properties of a spacetime described by the Einstein- 
matter system (II. ip . (II. 2p . (II. 3p . The Penrose singularity theorem men- 
tioned above is of this nature, and part of our arguments will also be pre- 
sented in this form. 

However, in order to verify such general conditions, in particular with 
respect to the existence of local or global solutions to the corresponding 
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initial value problem, a specific matter model must be chosen, and in the 
present paper this is a collisionless gas. All the particles in the gas are 
assumed to have the same rest mass, normalized to unity, and to move 
forward in time. Hence, their number density / is a non-negative function 
supported on the mass shell 



PM := [gapp'^p^ = -1, / > 0} 



a submanifold of the tangent bundle TM of the spacetime manifold M; p" 
are the canonical momenta corresponding to general coordinates x" = {t, x"") 
on M. We use coordinates {t, x"") with zero shift, and Latin indices run from 
1 to 3. On the mass shell PM the variable p^ becomes a function of the 
remaining variables {t,x"',p^): 

The number density / = f{t,x°',p^) satisfies a continuity equation, the so- 
called Vlasov equation, which says that / is constant along the geodesies of 
the spacetime metric, 

dtf + ^ d^af - 1 r^,pV d,af = 0, (1.4) 

where F^^ are the Christoffel symbols induced by the metric ga/s- The 
energy-momentum tensor is given by 



,1/2 dp dp dp^ 
PaPpfm ' , (1-5) 



where \g\ denotes the modulus of the determinant of the metric. The system 
p.ip . (jl.4p . (jl.Sp is the Einstein- Vlasov system in general coordinates. For 
an introduction to relativistic kinetic theory and the Einstein- Vlasov system 
we refer to [l] and [33] . 

If, for comparison, the matter is to be described as a perfect fluid with 
density IZ, four- velocity field C/", and pressure P, then the matter evolution 
equations are the Euler equations 

v^Ve^n + (7^ + p)v"Ua = o, 

(7^ + PWVc^Up + {gaf3 + UaUp)V''P = 0, 

where is the covariant derivative corresponding to the metric gafB- The 
energy-momentum tensor in this case is 

Tal3 = nUaUp + P{gat3 + UaUfs). 
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To close the Einstein-Euler system it has to be supplemented by an equation 
of state P = P{TZ). The choice P = yields the dust matter model referred 
to above. 

Due to the complexity of the field equations (jl.ip very little can be 
said about the questions at hand for these equations in their general form. 
Since on the other hand these questions are of considerable interest also 
in spacetimes satisfying simplifying symmetry assumptions, we from now 
on focus on asymptotically flat, spherically symmetric spacetimes and write 
down the metric 

in Schwarzschild coordinates. Here t G M is the time coordinate, r G [0, oo[ 
is the area radius, i.e., 47rr^ is the area of the orbit of the symmetry group 
S0(3) labeled by r, and the angles 6 £ [0, vr] and ip G [0, 27r] parameterize 
these orbits. Asymptotic flatness means that the metric quantities A and // 
have to satisfy the boundary conditions 

lim A(t,r) = lim /i(t,r) = 0. (1.6) 

For a metric of this form the 00, 11, and 01 components of the Einstein 
equations are found to be 

e-^^{2rXr - 1) + 1 = 8^Tr^e-'^^'Too, (1.7) 

e"2^(2r/x^ + 1) - 1 = Svrr^e-^^Tii, (1.8) 

Xt = 47TrToi, (1.9) 

where subscripts indicate partial derivatives. The 22 and 33 components are 
also nontrivial, but they are not needed for our analysis, and the remaining 
components vanish identically due to the symmetry assumption. 

Our aim is to find explicit conditions on the initial data such that the 
corresponding solutions of the spherically symmetric, asymptotically flat 
version of the system (jl.ip . (jl.2p . (jl.3p have the following property: There 
is an outgoing radial null geodesic 7"*" originating from r = tq > 0, i.e., 

^(s) = e{/^-A){s,7+(s))^ ^+(0) = ro, (1.10) 
ds 

such that the solution exists on the outer region 

D:={{t,r)G [0, cx)[V > 7+(i)}, (l-H) 
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and 7''' has the property that 

hm 7+(s) < oo. (1.12) 

This indicates that the matter distribution undergoes a gravitational col- 
lapse, and a black hole forms. In the case of Vlasov matter we obtain a 
more detailed picture which supports this interpretation: There exists an 
extremal, radially outgoing null geodesic 7* in the outer domain D such 
that lims^oo7*(s) = 2M where M is the ADM mass of the solution, and 
as t — > cxD the metric converges for r > 2M to the Schwarzschild metric 
representing a black hole of mass M. 

In the next section we state our main results for the Einstein- Vlasov 
system, where we specify classes of spherically symmetric initial data which 
lead to solutions showing the above behaviour. The Vlasov equation and the 
corresponding energy-momentum tensor components in the case of spherical 
symmetry are stated there. In Section [3] we give a general formulation of 
one of our results where no particular matter model is considered. The 
reason for this is that most steps in the proof of Theorem 12.21 below are of a 
general character and — besides the fact that for the Einstein- Vlasov system 
there is an existence theory for the initial value problem which guarantees 
the existence of solutions on D — the specific properties of Vlasov matter are 
used only in one key lemma. Hence it is natural to precisely single out the 
required conditions on the level of the macroscopic matter quantities. This 
clarifies the main mechanism in our method, and it may lead to applications 
of our method to other matter models. Using an additional feature of Vlasov 
matter we construct an alternative, and in some respects larger, class of 
initial data which ensure the formation of black holes, cf. Theorem 12.11 

The proofs of our results then proceed as follows. After stating some 
general auxiliary results in Section [J] we prove Theorem 13. H which is the 
general-matter version of Theorem 12.21 in Section El The latter result is 
then established in Section [6] by showing that Vlasov matter satisfies the 
required general conditions on the matter for a suitable class of initial data. 
Theorem 12. II is established in Section [7] together with Corallarv l2.3l on global 
existence in Schwarzschild coordinates. For all these results it is essential to 
make sure that in the outer region D all the matter moves inward. In the 
case of general matter this is a condition which we have to impose on the 
solution, whereas in the case of Vlasov matter we can specify conditions on 
the initial data such that this is true. In Section [8] we prove the convergence 
of our solution to a Schwarzschild black hole of the corresponding ADM 
mass in the case of Vlasov matter. 
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2 Main results for Vlasov matter 



In this section Eqns. (jl.6p - ()1.9p will be supplemented by the spherically 
symmetric version of the Vlasov equation together with expressions for the 
relevant components of the energy-momentum tensor so that a closed sys- 
tem is obtained, known as the spherically symmetric, asymptotically flat 
Einstein- Vlasov system. 

In order to exploit the symmetry it is useful to introduce non-canonical 
variables on momentum space and write / = f{t,r,w,L). For a detailed 
derivation of the corresponding equations we refer to [23]; here we just state 
the result. The Vlasov equation is 

dtf + e^'-^^drf - (^Xtw + e^'Vri? - d^f = 0- (2-1) 

where 

E = E{r, w, L) := Vl + 'W^ + L/rS = eV- 

The variables w g] — oo, oo[ and L G [0, oo[ can be thought of as the radial 
component of the momentum and the square of the angular momentum 
respectively. Notice that the latter is conserved along characteristics of the 
Vlasov equation. The matter quantities are given by 

p(t, r) = e-2^roo(t, r) = - / Ef{t, r, w, L) dL dw, (2.2) 
p{t,r) = e-^^Tu{t,r) = - / —f{t,r,w,L)dLdw, (2.3) 

i/ — ■ OO t/ 



j(t,r) = -e-(^+^)Toi(t,r) = ^ / w fit,r,w, L) dL dw. (2.4) 

^ J-oo Jo 



Notice that the quantities p, p, j appear on the right hand sides of the 
field equations (|1.7p ~ (|1.9p . and they are given in terms of / alone, which 
is the main reason for using the non-canonical variables w and L. The 
system (jl.6p - (|1.9p . (j2.ip ~ (j2.4p is the spherically symmetric Einstein- Vlasov 
system in Schwarzschild coordinates. As initial data we need to prescribe an 
initial distribution function /= f{r,w,L) > 0, which should be compactly 
supported in ]0, oo[x] — oo, oo[x]0, oo[, and such that 

l*r pr poo pQC ^ ^ 

/ 47r?7^p(r/) dry = 47rM / / Efi7],w,L)dLdwdr]<-. (2.5) 
Jo Jo J-oo Jo 2 

The origin r = is excluded from the support for technical reasons, but this 
could be avoided by using Cartesian coordinates. The condition ()2.5p implies 
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that the equations (|1.7p and (jl.Sp have solutions A and /x, cf. Section HI and 

o 

since / has compact support, a property which is inherited by f{t), the 
matter terms are well defined. If in addition / is C we say that the initial 
data is regular. As is shown in or [24], regular initial data launch a unique 
local solution for which all the derivatives which appear in the system exist 
classically. In Section [6] we discuss in more detail that this local solution 
extends to the whole outer region D defined in (jl.lip . 

To state our main results let < tq < ri be given, put M = ri/2 (this 
is going to be the ADM mass of the solution), and fix < Mout < M such 
that 

— < -. (2.6 

ro 9 

Remark. The value 8/9 is chosen for definiteness, and any number less 
than one would do, effecting the values of some of the constants below. 

Two different theorems will be stated below, corresponding to the fol- 
lowing two situations. 



(i) Let Ri > ri be such that 
or 

(ii) let Ri > ri be such that 



i^i-n<^, (2.7) 



Ri-n . Jl r^ ri - rp 



Ri [6 UkRiM' 8kRi 

where the (explicit) constant k > will be specified in Theorems 
and 13.11 below. 

Finally, we define 

Ro := ^(ri + i?i). 

Denote by p the energy density induced by the initial distribution function 
/. We require that all the matter in the outer region [ro,oo[ is initially 
located in the strip [Rq,Ri], with Mout being the corresponding fraction of 
the ADM mass M, i.e.. 



roo i-R\ 

/ ATTr'^p{r)dr = / Airr'^ p(r)dr = Mont- (2.9) 
J ro J Ro 
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Furthermore, the remaining fraction M — Mont should be initially located 
within the ball of area radius rg, i.e., 



/ 4Trr'^p{r)dr = M - Mont- (2.10) 
Jo 

Remark. The set up described above is quite similar to the set up in [9j 
for a scalar field. In [9] it is not required to have matter in an "inner" strip 
[0,ro], as is the case here in view of ()2.10p and the condition Mout < M. 
The reason why we need some matter in the region r < ro is to ensure that 
initially ingoing matter continues to be ingoing for all times, cf. Lemma |6. II 
below. If one only considers purely radially ingoing particles, i.e., with no 
angular momentum (which results in a non-smooth distribution function /), 
then we could allow for Mout = M. It is interesting to note that p = p holds 
for Vlasov matter, if the particles have no angular momentum and their rest 
mass is zero, which is the case for the scalar field considered in [9j. 

Now we are in the position to formulate our main results for Vlasov 
matter. Corresponding to Case (i) above, we prove 

Theorem 2.1 Let ro, ri, M, and Mout be given as above, and let Ri satisfy 



(2.1). Then there exists a set X\ of regular initial data for the spherically 

o 

symmetric Einstein- Vlasov system such that if f & Ii, then 112. 9\) and 112. 10\) 
hold, the corresponding solution exists on D, and 

/•OD 

lim 7^(s) < oo, lim / Airr'^ p{s , r) dr > 0, 



where 7"*" satisfies 111.10]) . 



By abuse of notation we denote by D both the outer region in spacetime 
defined by (jl.lip and the part of the mass shell with {t, r) G D. 

The next theorem addresses Case (ii) above, assuming the stronger condi- 
tion (12. Sp . This allows for a more straightforward proof, and the constraints 

o 

on the momentum variables of the initial distribution function / which are 
used to specify the set Ii will be slightly relaxed. Hence, the initial data set 
Ii does not contain X2 in Theorem 12.21 below, but it is larger in the sense 
that data in I2 are quite close to containing a trapped surface, which is not 
necessarily the case for data in Zi. The precise form of Ti and Z2 is specified 
in the proofs. 

Theorem 2.2 Let ro, ri,M, and Mout be given as above and let Ri satisfy 
\2. 8\) with K, = 6. Then there exists a set I2 of regular initial data for the 
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spherically symmetric Einstein- Vlasov system such that if f £ I2, then \2. 9^) 
and (2.10]) hold, the corresponding solution exists on D, and 

/•oo 

lim ■y~^{s) < 00, lim / 47rr^ p{s , r) dr > 0, 
where 7+ satisfies 111.10]) . 

The Einstein- Vlasov system has a wide variety of static, spherically symmet- 
ric solutions with finite ADM mass and finite radius, i.e., compact support 
of the matter, cf. |23l [271 130j . Particularly interesting examples of initial 
data for which our results apply are obtained if the matter for r < tq is 
represented by such a static solution, more precisely: 

Corollary 2.3 Let fs be a static solution of the spherically symmetric 
Einstein- Vlasov system with finite ADM mass Mg > and finite radius 
rg > 0. Define r^ = rg, let ri > tq be arbitrary, M = ri/2, and 
Moy^t = M — Ms', the latter quantity is positive. Then the initial data sets 

o 

Ti and T2 both contain data f which coincide with the given static solution 
for < r < rg . The corresponding solution f of the Einstein- Vlasov system 
exists for all r > 0, t > and coincides with the static solution fg for all 
r < 7+(t) and t > 0. 

We prove this result at the end of Section \7\ It represents a global 
existence result for the Einstein- Vlasov system in Schwarzschild time for 
data which are not small. 

In the next section we formulate a version of Theorem 12.21 for quite gen- 
eral matter models. One reason for this is that the main mechanism behind 
our method becomes very transparent by posing sufficient conditions on the 
macroscopic matter terms rather than conditions on the initial distribution 

o _ 

function / as we did in the theorems above. Theorem 12.21 will then be a 
consequence of this generalization, cf. Section [6l whereas Theorem 12.11 is 
established in Section [71 

In these proofs it turns out that for the classes of initial data that we 
specify we can obtain somewhat sharper asymptotic information on 7"'" and 
the mass in the outer region; see (15. Sh below. More importantly, we can 
establish the following additional information which shows that the solution 
evolves towards a Schwarzschild black hole of mass M. 

Theorem 2.4 In the situation of Theorem \2.1\ or Theorem \2.2\ the following 
holds: 
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(a) There exist constants a,(3 > depending only on the initial data set 
Ii or I2 respectively such that if 

t>0 andr>2M + ae~^* 

then f{t,r,-,-) = 0, i.e., we have vacuum, and the metric equals the 
Schwarzschild metric 

2M\ ,,,, . 2M\"^ , 9 , 2/ ,/,2 



ds^ = -[^- —j + - —j + ^ (^^ + ^'^'^ )' 

representing a black hole of mass M. 
(b) For allt>0 and 7+(i) <r<2M + ae'l^^ 



n{t, r) < In 



2M + ae-f^t 



so that in the outer region D, 

lim fi(t, r) = —00 for r < 2M, 



and the timelike lines r = c, where c G [0, 2M], are incomplete and 
their proper lengths are uniformly bounded by a constant depending on 
a, (3 and M. 



(c) Let 



r* := sup{r > ro | the radially outgoing null geodesic 7 with 
7(0) = r satisfies lim 7(3) < 00}, 



and let 7* be the radially outgoing null geodesic with 7* (0) = r* . Then 

lim 7*(s) = 2M, 

s— >oo 

and every radially outgoing null geodesic 7 with 7(0) > r* is future 
complete with linis—,oo"y{s) = 00. 



3 The result for general matter models 

In this section we specify the general assumptions on a matter model suffi- 
cient for our method to be applied. In order to keep the discussion consistent 
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with the Vlasov part of our arguments, and in view of the right hand sides 
of the field equations (jl.7p . (jl.Sp . (jl.Op . it is convenient to use the notation 



P ■= 



00, p:=e Til, j ■ 




(3.1) 



Firstly, we assume that the following two conditions are satisfied. 



• The dominant energy condition holds. 



(DEC) 



• The radial pressure p is non-negative. 



(NNP) 



The dominant energy condition (DEC) plays a central role in general rel- 
ativity and is the main criterion that a matter model should satisfy to be 
considered realistic. We refer to [l6j for its definition. The non-negative 
pressure condition (NNP) is restrictive in the sense that it rules out, for 
example, a Maxwell field as matter model. However, for most astrophys- 
ical models it is a standard assumption, with e.g. fluid models satisfying 
this condition. For the purpose of this paper we only need to focus on two 
consequences of these two criteria, cf. [16] and ^22j. The (DEC) condition 
implies, together with the (NNP) condition, that 



Furthermore, by (DEC) any geodesic {s,R{s)) of a material particle or a 
light ray satisfies 



The meaning of the latter condition is that locally the speed of energy flow 
is less than or equal to the speed of light. 

Let A, fi, p, p, j correspond to a solution of the spherically symmetric 
Einstein- matter equations ()1.6p - ()1.9p . (II. 2p . ()1.3p in Schwarzschild coordi- 
nates, launched by initial data from a class I. In order to investigate the 
global structure of the solutions it is necessary that they exist globally in 
an appropriate sense. In the situation at hand they need to exist on the 
outer region D defined in (|l.lip . In the spherical symmetric case the main 
obstruction for obtaining global solutions arises from the difficulties related 
to the centre of symmetry r = 0. For example, for a massless scalar field or 
a collisionless gas as matter model it has been shown that solutions remain 
regular away from r = for general initial data, cf. [121 [2], [28]. On the 
other hand, for dust a singularity of shell crossing type can also occur at 
some r > 0. Although in that case there are no true geometric spacetime 
singularities, such behaviour has to be ruled out in order not to interfere 




(3.2) 



dR{s) 



< g(/.-A)(s,R(s))^ 



(3.3) 



ds 
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with the analysis of the solution on D. This can be achieved by proper as- 
sumptions on the initial data, cf. [7j . In view of (j3.3p a possible break down 
of solutions at r = will have no influence on the outer domain D. Hence 
we formulate a third condition, concerning global existence of solutions in 
the outer domain, as follows. 

• For solutions launched by data from the set X, 7+ defined by (|1.1U|) 
exists on [0,oo[, and A, /i, p, p, j G C^{D). (GLO) 

The three conditions above are of a quite general nature. The fourth and 
final condition however, is tightly connected to our method of proof. 

• There exists a constant ci > such that p < —cij in D. (GCC) 

The acronym (GCC) stands for "gravitational collapse condition" , and this 
condition plays a crucial role for our method of proof. We emphasize that 
our main results show that for Vlasov matter there are initial data sets such 
that (GCC) holds. As a first consequence of (GCC) and (13. 2p . note that 
J < in D, i.e., the matter is ingoing for all times. In this respect our 
present results complement [4j, where purely outgoing matter is considered. 

Let us now assume that our matter model satisfies (DEC) and (NNP), 
and that there exists an initial data set 2 such that (GLO) and (GCC) hold 
as well. Then we have the following result, which should be viewed as a 
version of Theorem 12.21 for general matter. 

Theorem 3.1 Let tq, ri, M, and M^ut be given as above and let Ri satisfy 
\2. <g)j with K = 2ci. Assume that there exists an initial data setl^ C X such 
that ^2. 9(1 and ^2.10\) hold for all initial data in X3. Then for any solution 
launched by initial data in X3, 



where satisfies U.10\) . 

The detailed information on the gravitational collapse which for Vlasov mat- 
ter is provided in Theorem 12.41 is not available in the present situation, but 
the following still holds. 

Remark. In the situation of Theorem 13.11 



lim 7'^(s) < 00 



00 




lim p{t, r) 



00 for lim 7~^(s) < r < ri 



s— >oo 
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for some ri > linis^oo 7^(s)- If f* and 7* are defined as in Theorem 12.41 



and every radially outgoing null geodesic 7 with 7(0) > r* is future complete 
with lims^oo7('S) = 00. 

These assertions will be established in Section [8l Concerning the ques- 
tion which matter models besides Vlasov matter satisfy our conditions above 
we note the following: 

Remark. For a spherically symmetric perfect fluid with density TZ, pressure 
P = PijV), and radial velocity field tt, the (DEC) and (NNP) conditions and 
Eqn. (j3.2p respectively are satisfied provided that < PiJZ) < TZ, which 
restricts the equation of state. The (GCC) condition holds for example with 
ci = \/2 if —e^u > 1 on D. In the kinetic context of the Vlasov model 
we derive analogous estimates on the particle level from conditions on the 
initial data. 

4 Preliminaries 

In this section we collect some general facts concerning the spherically sym- 
metric Einstein-matter equations under the assumptions (DEC) and (NNP) 
that have been specified in the previous section. 

A quantity which plays an important role is the quasi-local mass m{t, r). 
Typically, the spherically symmetric Einstein-matter system is supple- 
mented by the requirement of a regular centre, i.e., A(t,0) = 0. Using 
this boundary condition the field equation (jl.7p implies that 



where the quasi-local mass would be given by m{t,r) := ^TTrf^ p{t,r]) dr]. 
Then m(t, 00) is a conserved quantity, the ADM mass. However, in the 
present context we want to investigate the system on the outer domain 
D, regardless of whether or not the solution remains regular in the region 
r < 7'^(t). Hence we do not use the usual boundary condition at r = 0. 
Instead, we assume that the ADM mass M > is given and redefine the 
quasi-local mass by 



then 



lim 7*(s) < 00 



s 



■00 




(4.2) 
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Then limj.^00 rn{t, r) = M, < m < M, and rrir = ATrr^p holds. Defining 
A by (j4.ip . (j3.ip shows that (jl.7p and the boundary condition in (|1.6p are 
satisfied. In addition, we need to modify (j2.5p to 

rn(r) < ^, rG]0,oo[, (4.3) 

a condition that once again will be included in the notion of regular initial 
data. 

By dLZl) and (fTHD . 

= [Attvp - ^)e2^, pr= + 47rrp)e2^. (4.4) 

In view of (|1.6p . p = fi + fl, where we define 

A(t,r):=- r^e2^(*'^)d^, (4.5) 

/•oo 

il{t,r):=- ATTr]p{t,r])e^^^^'''^ dr]. (4.6) 



Lemma 4.1 T/ie following assertions hold, 
(a) 2ji < ji — \ < ji < jl + \. 
(h) p + X<fi + X. 

(c) ifi - A)(t, r) = 2fi{t, r) + /J^ Anr, {p - p){t, r/) e^^^*' dr/. 

("d; /it(t,r) = /^°°47rj(i,?/)e(^+^)(*''')e2-^(*''') d??. In particular, if j < 0, 
i/ien a/so fit <0. 

Proof : In view of (|1.6p . 

[Attt] p{t, r/) - j '^'^ ^ ~ ] ^'^'^ ""^^ '^'^ " ^' 

and by (j3.2p the relation p — \ > 2fi follows. On the other hand, by ()4.ip . 
A > 0. Thus ft < leads to p — X<p<fL<fL + X, and part (a) is 
established. Part (b) follows from fi <0. As to (c), we observe that 

Attt] {p - p) dri = p, 

which gives the claim. By and ^M), {e^^^)^ = i^ie^^-'^)t = 

-A-K e^'^^e^^j. Hence (d) follows from (fi3|) . □ 
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Lemma 4.2 For r G [0, oc[ the following holds: 

47r7?(p + p)(t,77)e('^+^)(*''')e2^(*''^)dr/ = l-e('^+^)(*'^) < 1, 



Proof : It suffices to integrate 

5^(e/^+^) = ef'+^ifir + Xr) = e^+^47rr {p + p), 



= Airrpe'^+^e'^^, (4.7) 

observing that linir^oo A(^) ''') = lim^^oo r) = lim^^oo r) = 0. For 
Vlasov matter, the first relation has been used in [21 Lemma 1]. □ 

Next we consider outgoing and ingoing radial null geodesies 7^ and 7", 
respectively. 

Lemma 4.3 Let 7^ be the solutions to 



ds 

Then 



(s) = ±e('^-^)(''^^('*», 7+(0) = ro < n = 7-(0). 



(a) 7^ is strictly increasing, s m{s,'y~^{s)) is increasing, and the limits 
lims^oo 7'''(s) G]ro,oo] and lims^oom{s,^^{s)) G [m(0,ro),M] exist. 

(h) 7~ is strictly decreasing, s ^ m{s,"f~ {s)) is decreasing, and the limits 
lims^oo7~(s) G [0, ri[ and\\m.s^oom{s^^~ {s)) S [0, m(0, ri)] exist. 



(c) The relation 

d_ 

ds ' 

holds. In particular, if j < and p = j = along 7^, then also 
^(A + A)(s,7±(s))<0. 



^(A + X){s,j^is)) = (pt - 47rr e^^\j ^ p) 



(t,r)=(s,7±W) 



Proof : Differentiating (|4.ip w.r.t. t and using (jl.9p implies that mt 
— 47rr^e^~'*'j. Since p > j according to (j3.2p . this yields 

— m(s,7+(s)) = mt(s,7+(s)) + m,.(s, 7+(s))-^(s) 

= (_We^-^j + 4vrr2pe^-^)|(,^^)^(^_^^(^))>0. 
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Thus part (a) is obtained from m < M. Since p > —j, the proof of (b) is 
analogous to (a). As to (c), note that by definition of /i, (jl.7p . and (|1.9p . 

^(/i + A)(.,7^(s)) 



ds 



( - - "7 ^ X ' 

V * + f^^~^ + At + Xr- 



(t,r)=(s,7±(s)) 



/it-47rr e^+^(jTp) 

{t,r-)={s,7±(s)) 

as desired. The last claim follows from Lemma 14.1( d). 



□ 



5 Proof of Theorem 13.11 

In this section we use the hypotheses stated in Section [3] to prove Theo- 
rem 13.11 The proof is short and emphasizes that the crucial mechanism is 
captured in the (GCC) condition. Our main results which show in particu- 
lar that the (GCC) condition holds for Vlasov matter are established in the 
next sections. 

Consider the out- and ingoing null geodesies 7"^ and 7" defined in 
Lemma 14.31 The claims follow if we can show that these geodesies never 
intersect. By continuity and monotonicity there exists T s]0, oo] such that 

ro<7+W <7"W <n, *G[0,T[; (5.1) 

it will be shown that actually T = oo holds. In view of (j2.9p we have initially 
that p = p = j = for r > Ri. The (GCC) condition implies that j < in 
D, meaning that the flow of matter is ingoing. Therefore 

p = p = j = and m = M for {t,r) e [0,T[x[Ri,cc[. (5.2) 

By Lemma|421 (IS^D, the (GCC) condition, and LemmaSHd) for s G [0, T[ 
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and r G [7^(s), oo[, 

j-OO 

l_e(A'+^)(«'^) = / 47rr?(p + p)(s,r?)e(^+^)(^''')e2^('^''')(ir/ 

J r 

/•oo 

<2ci / A^Tr]\j{s,7])\e^^'+^^^''''^e'^^^'^''Ur] 

J r 

/oo 

= -2cii2i/it(s,r), 
since j{s,ri) / implies 77 < Thus 

^*^''"^-"2^(^"'^''^'^^''''0- ^^-^^ 
This in turn imphes that 



A(t,7^(t))-/i(0,7^(0)) 
d 
ds 



fi{s,j^{s))ds 











- 2cii?i+yo V2ciiii+ 7±(^)2 ^ ^ ^ ^ 

Now for any r G [rg, ri] and t G [0, T[ it follows from /i,. > and (|4.ip that 



Using M = ri/2 we get 



1 — ri 

Kt,Ri) = o log 



2 ° V iii 
so that for r £ [ro,ri], 



^Kt,r)^^Kt,R^) = J^l^, (5.6) 

R^ 
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By ()3.3p and the properties of the initial matter distribution there is vacuum 
in the region ^^{t) <r<^~[t). Hence m{t,r) = M — Mout and (j2.6p imply 
that 



< 3 



(5.7) 



V^l - 2(M - Mout)/ro 

for j~^{t) < r < 7^(t). From Lemma l4.ir b) and ()2.8p . recalhng k = 2ci, we 
obtain in particular that 



1 



2' 8cii?iM 



Thus ([531) yields 

A(t,7^(i))-/i(0,7^(0)) < 



< 



< 



2cii?i 
f l-d 

1 

4cii2i 



+ d 







1 M 
2cii?i + 72 



SciRi 

Hence Lemma l4.1f a) leads to the estimate 

ft 



, tG[o,r[. 



|7^(t)- 7^(0)1 



< e 



g(M~-A)(s,7±{.))^, 
A(0,7=^(0)) 







< 







Ri-n 



^1 



where we used (|5.6p in the last inequality. By the third condition in (j2.8p . 



i?i - ri ^ ri - ro 



16cii?i 



so that 



|^±(t)_^±(0)| <!1_I0, tG[0,T[. 



Since 7~(0) — 7^(0) = n — tq, this implies that 7^(T) — 7^(T) > 0. Hence, 
if we choose T in (|5.ip to be maximal, then T = 00, i.e., 7"^ and 7" do never 
intersect. This completes the proof of Theorem 13.11 □ 
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Remark. In the above proof we have obtained the more exphcit information 
that 

hm 7+(s) < ^^^i^, m(s,7+(s)) = M-Mout, s>0, (5.8) 

the latter since all the matter originally to the right of 7~(s) > 7^(s) 
necessarily stays there. 



6 Proof of Theorem [22] 

We first check that the (DEC), (NNP), and (GLO) conditions hold for Vlasov 
matter. Then we show that there exists a class of initial data such that the 
corresponding solutions satisfy the (GCC) condition with ci = 3. Hence 
Theorem 12.21 will follow from Theorem 13. li 

The characteristic system associated to the Vlasov equation p. II) is 



^^^i^-X)is,R)W^ (6.1) 
ds E 

= -A,(,,fi)H'-el''-^)l'-'«^,(s,fi)£ + el'-^)C-'"-^, (6.2) 

as 

If s 1-^ {R, W, L){s) is a solution with data {R, W, L){0) = (r, w, L), then 
f{s,R{s),W{s),L) = f{r,w,L) 

o 

is constant in s. Hence {R{s),W{s), L) G supp/(s) iff {r,w,L) G supp/. 
Such characteristics will be addressed as characteristics in supp/. 

Direct inspection of the definition in (j2.3p shows that (NNP) holds for 
Vlasov matter. It is moreover well-known that the (DEC) condition is sat- 
isfied for Vlasov matter; see [H Sec. 1.4]. Alternatively, we can check (j3.2p 
and (j3.3p directly. The latter follows from (j6.ip above, whereas the former 
is a consequence of the expressions for the matter terms given in (j2.2p , (j2.3p , 
and ([23]). 

To see that the (GLO) condition holds for any regular initial data set 
we argue as follows. First of all, a regular initial data launches a local-in- 
time solution on some time interval [0, T[, and the corresponding theorems 
in [25] or |24j also give a condition under which this local solution can be 
extended to a global one. In order to see that the local solution can always be 
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extended to the whole outer domain D we first observe that the sphericahy 
symmetric Einstein- Vlasov system on D, with (j4.ip and (j4.2p replacing the 
usual boundary condition of a regular centre and with (jl.lOp included, has 
again a well-posed initial value problem for regular data supported in ]ro, oo[. 
This can be shown in the same way as for the system on the whole space, the 
essential point being that no characteristic of the Vlasov equation can enter 
region D at the boundary r = 'y~^{t). To the local solution on D we can now 
apply the arguments from [28] and conclude that the solution exists on all 
of D. This is possible due to the fact that the estimates in [28] address a 
situation where matter is bounded away from the centre or is controlled in a 
neighborhood of the centre so that these estimates can be applied on D. We 
emphasize that for our present analysis only the behaviour of the solution 
on D plays a role. We have chosen to present our results in the form that 
we have Vlasov matter also inside r < ^~^{t), and this part of the solution 
may or may not break down, but this is irrelevant for our arguments. 

Hence it remains to show that the (GCC) condition holds. To this end 
we let < ro < n < Ri, Rq = (ri + Ri)/2, and M = ri/2. For a parameter 

o 

W- < to be specified below and regular data / with ADM mass M we 
formulate the following 

o 

General support condition: For all (r, w, L) G supp / the following holds: 

r G]0,ro] U [i?o,i?i], 

and if r G [Rq, Ri] then 

w<W- 

and also 

0<L< — -mTr/) -I- jyrrifr/), G [ro, (6.4) 



We use the notation mwhen p = p m. (j4.2p . Furthermore, we abbreviate 



T = T{n,R,):=J^—^. (6.5) 

The following lemma shows that if the support condition holds, then the 
particles in the outer domain D keep moving inward in a controlled way. 

o 

Lemma 6.1 Let f be regular and satisfy the general support condition for 
some W- < 0. Then for all {r,w,L) £ supp/(t) such that {t,r) G D, 

w < T{ri,Ri)W-. 

In particular, j < on D. 
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Proof : Let [0, T[ denote the maximal time interval such that for t < T 

< for {r,w,L) £ supp/(t) with {t,r) £ D. (6-6) 
Since W- < 0, T > by continuity. By the definition of j, 

j{t, r) < for (t, r) £DT:=Dn ([0, r[x [0, oo[). 
Let {R,W, L){s) be a characteristic in supp/. Then 



(6.7) 



e~^{WXt + WXr I 



ds 



ds J 



inR 



e''{2WEj - W^p - E'^p) + ( 1 



2m \ L 



m /w 



+ e>'^ —-E 



i?2 V E 

47r2 
R 



oo />oo 



oo Jo 



'E IE 

— W — \ — W 

E \ E 



R J R^E 



f dL dw 



m (l + L/R^ 2L \ 



i?2 



E 



+ 



R^E 



L 



R^E 



where E = E{R, W, L) and E = E{R, w, L). Therefore 



m (l + L/R^ 2L \ 



E 



R^EJ 



L 



R^E' 



Differentiating (|4.ip w.r.t. t and using (jl.9p leads to mj = —A7rr'^e^~'^j, 
which by (j6.7p is non-negative on D^^. It follows that m{s,r) > m(0, r) = 
m{r). Thus as long as the characteristic remains in Dt-, 



(l + L/R^ 2L \ 



R^E 



E 
3L 



+ 



L 



R^E 



L - — m{R)-Rm{R) 
R 



Now R{{)) G [i?o,/?i] and < by 1^ and ([62]) yields Ri > R{0) > 
-R(s) > 7'''('S) ^ '"0- Hence condition (16. 4p implies that, as long as the 
characteristic remains in Dt, £{e~^W) < 0, so that 

W{s) < e^('^'^(^))--^(0'-f^(o)) W-. 
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But A > 0, so < leads to 

W(s)<( mill e-^(°''')W„. 

\re[Ro,Ri] ) 

In view of (jlTT]) . 



and recalling (j6.5p it follows that 

< r(ri,i?i)iy„ < 

as long as the characteristic remains in Dq^. By the maximality of T in (j6.6p . 
T = oo, and the proof is complete. □ 

In order to specify the initial data set I2, let ro, ri, M, and Mout be 
given as in Section [2] and let Ri be such that (j2.8p holds for k = 6. We 
require that < satisfies the estimate 

r{ri,Ri)\W^\>l. (6.8) 

Then 

/I / is regular, satisfies (|2.9p . (j2.10p . the general support condition, 
and for (r, w, L) G supp/ with r G [i?0; -^i]; ^/L/tq < T 

(6.9) 

Consider now a solution / launched by initial data from this set. Condition 
(j6.8p and Lemma |6 . II implv that 

\w\>T{ri,Ri)\W^\>l on supp/nZ?, (6.10) 

and since L is conserved along characteristics, (16. 9p leads to Vl/t < 
^/L/ro < \w\ for all particles in supp / D D. Hence the definition (|2.2p 
of p implies that on D, 

p{t,r)<— / fdLdw + — / |u;|/dLdu; 
^ J-00 Jo T J -00 Jo 

+ ^ / / v^/rfdLdw 

T J -00 Jo 

^ f'OO /'OO 

<3-/ / |u;|/dLd«; = 3|i(t,r)|. (6.11) 
^ J -00 Jo 
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Accordingly, X2 satisfies the (GCC) condition with c\ = 3, and Theorem 
follows from Theorem 13. 1[ □ 

We briefly show that the set I2 is far from empty. Therefore fix < ro < 
n < Ro < Ri, M = ri/2, and < Mout < M such that Rq = {n + i?i)/2, 
(|2.6p . and (j2.8p are satisfied. Let < /i G C'^ have r-support in [ro — (5, ro] 
for some < 5 < ro/9, and let < /2 G have r-support in [Rq, Ri\. Fix 
the compact tn-support of /2 in ] — cxd, WJ\ with W- < such that ([6? 
holds, and fix its L-support in [0, L2] so that 



and 



L < (M-Mout)(Y + ^)' Le[0,L2], r?G[ro,fli 



<T{n,Ri)\W. 



Now take / = Afi + -B/2, where A > and B > are chosen such that 
(fTOl) and (pH)]) are satisfied. Note that rh{r]) > M - Mout for r/ G [ro, Ri], 
whence ()6.4p holds as well; thus the general support condition if verified. It 
remains to check (14. 3j) . If r G]0, ro — (5], then rn(r) =0. If r G [ro — S,Ro], 
then rn(r) < M — Mout yields in view of (|2.6p . 



If r G [iio) 00 [) then 



2m ^ 2(M - Mout) ^ ^ 
r ~ ro — 5 



2m 2M 
r i?o 

o 

since 2M = ri < Ro- Hence / is regular and has all the properties that are 
required in the definition of I2 ■ 

Remark. The set I2 has "non-empty interior" , in the sense that sufficiently 
small perturbations of initial data in the "interior" of this set belong to I2 as 
well, provided that the support is changed very little and M is left invariant. 
This is due to the fact that the various parameters entering into the definition 
of I2 are defined in terms of inequalities and hence can be varied. 



7 Proof of Theorem 12.11 

The set up is closely related to the set up in the proof of Theorem 12.21 As 
we saw above, the (DEC), (NNP), and (GLO) conditions are satisfied for 
Vlasov matter, and we will again construct an initial data set such that the 
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(GCC) condition holds with ci = 3. However, since this result relies on 
condition (j2.7p instead of (j2.8p . we cannot simply invoke Theorem 13.11 after 
the (GCC) condition has been verified; instead an additional step needs to 
be added to the proof. For this new argument a slightly stronger condition 

o 

on the momentum variable w needs to be imposed on supp /. We now require 
that < satisfies 

r(ri,i?i)2|W„|2>i^, (7.1) 
a 

where 

ro ri-ro 



d := min 
Then 



2' l2Ri 300i?i 



{o o 
/I / is regular, satisfies (|2.9p . (j2.10p . the general support condition, 

and for (r, L) G supp / with r £ [Rq, Ri], \/L/ro < l.| (7.2) 

The same construction as at the end of the previous section shows that this 
set is not empty, and the same remark as at the end of the previous section 
applies. 

Let / be a solution launched by initial data from Xi. It is clear from 
these conditions that Lemma [6.11 applies . and since 10/ d > 1, it follows that 
()6.10p holds as well. Thus the argument leading to p < 3|j| on D in the 
proof of Theorem 12.21 applies again. Hence, the (GCC) condition is satisfied 
with ci = 3. 

Consider the expression 



^ POO i-oo . ^2 , 

p{s,r) - p{s,r) = / - —j f{s,r,w,L)dL 



dw. 



-oo Jo 

Since E"^ > up' > T'^{ri,Ri)Wl by Lemma WT\ we get for r G [7+(s),i2i] 
from ^fLjr^^ < 1, 

it;2l„„l/L\2 2 
E-^ = ^^E^-^^-} = ^{^ + ^)<^<^J=:coE, (7.3) 

SO that 

p{s, r) - p{s, r) < cop(s, r). (7.4) 

After this preparation, we again show that the out- and ingoing null 
geodesies 7"'" and 7" do not intersect. We choose T g]0,oo[ such that (jS.ip 
holds. In this case we cannot rely on the smallness of as in the proof of 
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Theorem 13.11 so we need to control the evolution also when is not small. 
For this part the estimate (j7.4p is essential. 
We fix G [0, T[ by requiring that 

g(/i+A)(s,7±(s)) > for s G [0,tt[, e{A+A)(s,7±W) < d for s e [t^,r[. 

First we note that is well-defined, since by Lemma H^ c). 

Ag{A+A)(s,7±W) <Q_ 
ds 

Step 1: Consider s G [0, t^]; if = 0, then this step is omitted. For 
since {fi + X)r = ^irrpe^^ > by (jiT]) . Hence Lemma HU c) and fTil) yield 

j-OO 

{^i - A)(s,7^(5)) = 2A(5,7^(s)) + / 47rr? (/9 -p)(s,r/) e^^^^'") drj 

J7±(s) 

< 2/i(s,7±(s)) + ^ / 47rr?p(s,?7)e('^+^)(^''')e2^(^''')dr? 

" ^7±(s) 

<2A(.,7±(s)) + ^, 

where for the last estimate Lemma 14.21 has been used. 

Now we make the following observation: There is at least one charac- 
teristic {R,W , L){s) with ^(0) G [i?o,-Ri], which does not leave the strip 
[ri,i?i] during the finite time interval [0,T]. In fact, if at time t = T all 
characteristics had left the strip [ri,i?i] (and thus had entered the region 
r < ri), then m{T,ri) = M. From (j4.ip and 2M = ri it would follow that 
A(r, ri) = oo. However, this contradicts the (GLO) condition which holds 
for Vlasov matter. 

Since 7^(s) < ri < R{s), and since fir > 0, we thus obtain in view of 
Lemma |4.1 f a) that 

(/X - A)(s,7±(s)) < 2/i(.,7±(5)) + ^ < 2A(s, Ris)) + ^ 
<(/x-A)(s,i?(s)) + ^, sG[0,tt]. 

Next note that \W\ > 1 by (j6.10p . and hence due to (j6.ip and observing 
i?2 > ^2 > L, 

' ^ ~ V2 + Ty2 - 2 
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Therefore for all t G [0, t J] the estimate 



|7^(t)- 7^(0)1 



^g(/.-A)(s,7±(s)) 



< e d 



<-2e^ / R{s) ds = 2e^ iR{0) - Rit)) 



£0 



< 2e-{Ri - n) 



(7.6) 



is obtained. 

Step 2: Let t S [t^,T[; if = T, then this step is omitted. The arguments 
here are basically the ones presented in Section [5j The computation leading 
to ()5.4p is almost identical, and 

for ci = 3. By LemmaSUb), e(^+^)(«' < e(/^+^)(^' < d. Next we 
use the facts that m/r < 1/2, 7^(s) > tq, and the definition of d to obtain 
the estimate 

r-* / 1 In 



< ^ + 



< 



< 



2ciRi ' 




/ 1 - d 


-d± 

2ro 


V2cii?i 


( ' 


1 

-d 

2ro 





2rn 



it - tt) 
it - tt) 
{t-tf), tG[tt,T[. 



Hence by Lemma 14.1( a). 

l7^W-7^(tt)l = 



< 



< g(A+A)(tJ,7±(t?)) / e~^^ < 8cii?id. (7.8) 



Adding the contributions (|7.6p from Step 1 and (|7.8p from Step 2, the final 
estimate 

|7±(t) - 7^(0)1 < 2e"o/'^(i?i - ri) + Sciiiid 
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is obtained for all t G [0, T[. From (|7.3p and (j7.ip we have co/d < 1/5. The 
third condition on d together with (j2.7p thus imply that 

|7±(t) -7^(0)1 < 



As in the proof of Theorem lS.ll we conclude that 7^ and 7 do not intersect, 
completing the proof of Theorem 12.11 □ 

Remarks, (a) The sharper estimates stated in ()5.8p clearly hold also in 
this case. 

(b) The solution must necessarily enter the regime of Step 2, more precisely, 

hm e(^+^)(^'^*W) =0 

for both null geodesies. Otherwise, the monotonicity implied by Eqn. (|7.5p 
yields a positive constant c > such that e^'^^'^^^''''^ '^^-'^ > c for all time, 
and hence. 

Since no matter can cross the two null geodesies, 

/oo 
47r7?(2p - p)e^^dr] + 2/i(s, r) 



>2A(s,r) = -2 / 
r — 2rn(ro) 



^?^o) 1 
ry2 1 — 2rn(ro)/r/ 



In- 



fer r = 7^(s). If we insert this into the estimate for 7^ it follows that this 
quantity is bounded from below by a positive constant which contradicts 
the finite limits of 7^(s) as s — > 00. 

It remains to prove Cor. 12. 31 
Proof of Corollary 12.31 : Let fs be a static solution. By [3j, 2ms{r)/r < 
8/9 for r > where nig is the local ADM mass induced by fs- In particular, 
Ms < rs/2 < ri/2 = M, and (j2.6p holds. As described above we can now 

o 

specify the matter distribution for r > rg, and we obtain initial data / in Xi 
or in I2 which coincide with the given static solution for < r < tq . 

Since no matter travels from the outer domain D to the inner one where 
f ^ the only way the matter in the outer domain can affect the static 

solution is through the metric. Consider the time-independent version of the 
Vlasov equation (j2.ip . Dropping all the time derivatives we see that in the 
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remaining equation the factor e can be canceled. Therefore, the static 
Einstein- Vlasov system is formulated in terms the quantities /, A, and fir, 
but not fi itself. By ()4.ip and (|4.4p . A and /i^ are on r < j'^it) not affected 
by the matter in the outer domain D, and therefore / = fs,\,fir remain 
time-independent for r < 7+(t). □ 

Notice that the metric coefficient fi of course does change on the interior 
region, cf. Thm 12.41 fb). 



8 Proof of Theorem [231 

As a first step we estimate /x— A from below for r > 2M, using Lemma [4. 11 (a): 
(f. - \)(t. r) > 2A(t, r) = -2 H !!%^e2^('-')<i. 



7] {f] — 2m{t,r])) J J. r] (r] — 2M) 

r-2M 
= In , r > 2M. 

r 

Now consider any characteristic in the matter support, and let R{t) denote 
its radial coordinate. Then by Lemma |6. II and as long as R{t) > 2M, 

ds E - - R ' 

for initial data from the set Xi respectively I2 one can take C := 



T\W-\/ y 2 + T'^W'^ respectively C := l/\/3. Integrating this differential 
inequality we find that as long as R{t) > 2M the estimate 



-Ct > r^'^ ^ dr = Rit) - R(0) + 2M In ^f^l 

- A(o) r-2M ^' ^ ' R{0) - 2M 

>2M-i?, + 2Mln.^(^)-^^ 



R{0) - 2M 
holds, and hence 

R{t) <2M + {Ri - 2M)e2^(^i-2^-^*\ 

which proves the support estimate in part (a). Since all the matter, which 
has ADM mass M, is contained in the region where r < 2M + ae~^* =: 
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cr(t), the assertion on the metric follows. Moreover, for any r < a{t) the 
monotonicity of // with respect to r implies that 

r) < ^(t, ait)) = flit, ait)) = In [-^^^ J , 

which is the first assertion of part (b). The second follows immediately since 
the integral dt [g the proper length of a coordinate line of constant 

r, 9, and (p in the outer region D. This completes the proof of part (b). 

As to (c) we first observe that any radially outgoing null geodesic which 
enters the region r > 2M escapes to r = oo and is future complete, since by 
part (a) the metric on r > 2M + e where e > is arbitrary eventually equals 
the Schwarzschild one for which the asserted properties of the geodesies hold. 
Now consider the extremal geodesic 7*. If there existed some time t > 
such that J* it) > 2M, then by continuous dependence on the initial data 
the same would be true for all radially outgoing null geodesies with 7(0) 
sufficiently close to but less than r* . Hence such geodesies would escape to 
r = 00 in contradiction to the definition of r* . This shows that the extremal, 
radially outgoing null geodesic 7* has the property that limj^oo 7*(^) ^ 2M. 

It remains to show that the limit above cannot be strictly less than 
2M. To this end we consider a radially outgoing null geodesic as long as 
-fit) < ait) = 2M + ae"^*. Then 



ds \ ai-s) J ' 



and hence for any Q <tQ <t and as long as 7(t) < cr(t), 

l(t)<iito) + Ce-f'^l\ 

where the constant C > again depends only on the initial data set. Now 
assume that R* := hmt^oo7*(t) < 2M, choose to > such that R* + 
(7g-/3to/2 ^ 2M, and consider the radially outgoing null geodesic 7** with 
7** (to) = R*- Then by construction, 7**(t) < 2M < ait) for ah t > to, and 
since 7** (to) = R* > 7* (to) it follows that 7**(0) > 7*(0) = r*. Hence 7** 
is a radially outgoing null geodesic which at time t = starts to the right of 
r* and does not escape to r = 00. This is in contradiction to the definition 
ofr*. □ 

We conclude this section by proving the remark after Theorem 13. 11 Un- 
der our general matter conditions the matter is ingoing in the region D, in 
particular, the matter is for all time restricted to the region where r < Ri. 
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Hence for r > Ri the metric is again equal to the Schwarzschild one with 
mass M, and if we replace 2M by Ri in the above argument for part (c) 
we obtain the assertions on 7* in the general matter context. As to the 
divergence of we observe that 

^/i(s,7"(s)) = /it(s,7"(s)) +/ir(s,7"(s))^(s) < 

so that the limit /ioo := limg^oo /"(s, 7~(s)) exists. The fact that 7~(s) > 
is decreasing with 

= g{/^-^){*.7"(«)) > g2/i(s,7"(s)) > e^Aoo^ 

implies that /ioo = — 00. Since ^ < jl we conclude that 

lim ^(s,7~(s)) = —00, 

s— >oo 

from which the assertion follows by the monotonicity of fi with respect to r. 
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